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PHOTOPOLARIMETRY OF SCATTERING SURFACES AND 
THEIR INTERPRETATION BY COMPUTER MODEL 

SUMMARY 


The cc»iipuv>,'r model of Wolff (1975) of a rough planetary surface 
is simplified and revised. Close adherence to the actual geometry of 
a pitted surface and the inclusion of a new function for diffuse light 
has resulted in a quantitative model comparable to observations by 
planetary satellites and asteroids. 

A new function is derived to describe diffuse light emitted 
from a particulate surface. The function is in terms of the indices 
of refraction of the surface- material , particle size, and viewing 
angles. 

Computer-generated plots are made that describe the observable 
and theoretical light components for the Moon, Mercury, Mars and a 
spectrum of asteroids. Other plots describe the effects of changing 
surface material properties. 

Mathematical results are generated to relate the parameters of the 
negative polarization branch to the properties of surface pitting. 

The range of validity and accuracy of the various Albedo-polarization 
rules is discussed in terms of the computer surface model, and the CSM 
taxonomy is matched with the predictions of the Model. 

An explanation is offered for the polarization of the rings of 
Saturn, and the average diameter of ring objects is found to be 30 to 
40 centimeters. 
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light, 2) doiily-reflected light, 3 ) diffusely reflected light, 4) effects 
of shadow that creates the opposition effect, and 5) effects of snade on 
dodDle- reflections that create the negative branch. 

Despite the model's success in dealing with polarization, it suffered 
in demaiding long computer routines to examine even minor parameter 
changes, and was not subject to mathematical analysis to provide further 
insights . 

In this article, an Improved model is created in vdiich the detailed 
mechanists of reflection in a surface are exanined so as to provide the 
model with a closer relationship to the physical structure of the surface. 
Also, the long computer routines are replaced by a simpler model of the 
double-reflection process using four representative light paths to 
represent the doiJsles processes in each of four quadrants. The ’esult is a 
simpler, more tractable model which can be calculated easily and which 
provides additional insight to the reflection processes taking place in a 
rough surface. The model provides both polar imetric and photometric 
equations as a function of the phase angle and physical surface properties. 
The Model 

It is assuned that the solar-system environment contains a meteroid 
flux Wiich more or less uniformly impacts the surface of airless planets. 
The resulting surfaces are not assuned ideitical but the impaction process 
produces similar pits and jagged particles ranging fhom micron-sized dust 
to larger particles which are in turn covered with pits and particles. All 
variations of this type of surface are sufficiently similar to justify a 
model of light reflection in which the dominant processes are: 

1 ) single reflection 

2) double reflection in pits and between psarticles 
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3) diffuse, non- polarizing reflection, and, 

^0 the production of shade and shadow 

The identification of variations fhom the model in measurements of 
solar system bodies is an objective of the analysis. 

Diffuse reflection is defined as the sun of all processes except 
single and doiible reflection, and includes refraction, higher-order 
reflection, and particle scattering. It is assuned that the diffuse light 
is uipolarized because processes requiring three or more events in a random 
Jagged surface will almost certainly produce randonized light. 

The Unique Properties of Pits 

Despite the apparent disorder of a rough particulate surface, it 
turns out that there exist mathematically describeable properties of such 
surfaces, which are important to the photopolarimetric measurements of the 
reflected light. Unfortunately, perhaps because scientists are disciplined 
orderly thinkers, we have acquired an antipathy to the apparent disorder of 
a dirt surface. Reflecting this attitude, it has been ccmmonplace to 
sudbstltute for a surface, an array of s^eres or a cloud of random 
particles which can be attacked by neat, known mathematical methods. It is 
argued here that this has led to some erroneous conclusions, and, if one 
wants to have a reliable model and a descriptive mathematical formalism, 
then it is necessary to directly exanlne the geometric properties of the 
real surface itself. 

The random pits, holes and interstices of an impacted surface have 
the geometric property of partically trapping light rays in a manner that 
suppresses the probability of single-reflection while enhancing the 
probability of doiiJle- reflection. Consider the fate of a ray vhich enters 
a pit but has not yet headed outwards. The probability of getting out with 
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only a single reflection decreases the farther in the ray goes. The 
probability of a single- reflected ray getting out is equal to the solid 
angle occupied by ”clear sky” (as the ray sees things from down in the 
hole) divided by 4TT. That is 

Singles probability * P-| a (angular apeture)/4 7T 
Note that we are not talking about the probability of optical reflection 
which is given by the Fresnel equations and must also be a factor, but 
about the geometric probability of "getting out" or happening as a single 
reflection, as if the rays were billard balls. In that follows this will 
always be called the geometric probability (of the surface structure) . 

Note that if . singles geometric probability and ?2 is for doubles, etc. 
then 

+ P 3 .. =1 

We can approximate the singles geometric probability by assuning that the 
diameter of pits is about equal co the depth of pits. Then, we have that 
P-j sTTd^/^hrdS* 4 r 1/16 which is very much snaller than that obtained fhom 
cloud models. 

It is the character of a pit to retirn a ray back into its sourqe 
hemisphere with two or more reflections. Accordingly, the geometric 
probability of dotble reflection and higher order reflections are enhanced 
over cloud models or flat surfaces. Although pits come in all shapes, the 
concept is illustrated at opposition in Figure 1 using hemispherical pits. 
TWO reflections will always occur in Figure 1A if a ray enters inside of 
the center circle, so the geometric probability is proportional to the 
circle area. If a ray enters the annulus outside, three or more 
reflections will occur with a probability proportional to the annulus area. 

The relative probabilities are approximated knowing that the circle 
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FIGURE In This illustrates an optical 
property of pits. At opposition a hesd- 
spherical pit creates about four tines as 
many double reflections as singles. All 
higher multiples ire about equal to the 
amount of doubles, 

A more realistic surface behaves 
similarly to the hemispherical surface, but 
depends on the extent of the surface ruggedness. 
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P)) SINGLES B)FOR;VJftRO'OOUBUES c)BBC«WARD DOUBLES 

FIGURE 2. The angular distribution of the reflected light depends on the 
average vddth to depth ratio of the surface pits. The ang>alar distribution 
also depends on the type and direction of the reflection. Forward-going 
double reflections are most constricted in entrance and exit from a pit. 
Backward-going doubles are least constricted, with sideways doubles and 
singles being intermediate, ^ilD is a computer parameter proportional to 
the width to depth ratio. 
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radii are the sides of a 45° - 90° triangle. So, using the fact that the 
sun of all probabilities equal unity, we obtain 

P, 1/8 

Pg ^ 4/8 

Pj + P^ ... ^3/8 

we have arrived at the important result that the geometric probability of 
double reflection is about 4 times larger than single reflection. This is 
quite opposite to the Impression gained from using randan cloud models that 
ignore the special character of a pitted surface. Of course, we have not 
yet counted the single reflections that strike the external surface without 
entering a pit. If such external areas were assuned equal to the space 
between touching round pits, then the P^ geometric probability beccmes 
about 1/4 instead of 1/8. P 2 , P 2 ®tc. are decreased jroportionately. 

Relatively simple observations of the rough airless planets allows 
one to deduce that the proportion of singly-reflected rays is very small in 
the light emitted from them. Indeed, a convincing experiment can be done 
with two sheets of paper - one of them glossy and one of them rough. Look 
at the reflection of an electric lanp from the glossy paper, at all 
possible phase angles. One can easily see the large increase of 
reflectivity predicted by Fresnel's equations at phase angles of 160-180 
degrees. Note that this phenomenon identifie s single reflection. Now, 
repeat the experiment using the ro’.igh paper; there is no increase seen at 
large phase angles. It is concluded that single-reflections from the rough 
paper, and the rough planets, are highly attenuated - contrary to 
conclusions based on random cloud models. 


On the Effects of Fine Dust 


Electron microscope photographs of luiar surface particles by Bowell, 
DollfUs, and Geake (1972) showed the particles to be heavily coated with 
tiny particles of wavelength- si zed dimensions (1 micron). This dust has a 
major effect on the optical properties of the surface. 

Because of the small size of the dust particles and this top- layer 
location, the dust is able to reflect much more light than its proportion 
of mass in the total surface layer would indicate, especially at 
opposition. The opposition effect is a phenomenon due primarily to 

shadowing by the large particles on a surface (particles which are opaque 
to the light being observed). VIhen not in shadow, the dust particles 
adtering to the bottoms of pits enhance the reflectivity. If the dust 
particles are not opaque, which is to be expected if they are small, then a 
large fhaction of the light refhacted into them is reversed by internal 
reflections and emitted to the observer. The bottoms of pits and other 
areas are then much brighter than if the dust were not there. 

The presence of fine dust might be detected by measurement of the 
amplitude of the opposition peak. Knowledge of fine dus.t presence may be. a 
useful clue to the history and conposition of the body observed. 

Calculations of the computer model for the Moon and ccraparison with 
measured properties of linar rocks lead to the deduction that the Moon must 
have considerable fine dust. Otherwise, the reported average size (about 
20 microns) of particles is inccmpatible with the lack of specular 
reflection fhom the Moon. 

Because dust, if present, is the first reflecting layer , it is of 
prime interest. Especially so since the attention of research activities 
has tended to overlook the tiny volune of dust in favor of the more 


interesting large particles and rocks. It has become very easy to assume 
that the optical surface consists of the particles and rocks rather than 
the inconspicuous dust. 

The computer surface model can accommodate one surface layer which 
can be particles of any size (an average size is assumed). In the model, 
the uppermost layer gives rise to the diffuse component of light which is a 
function of particle size, indices of refraction, wavelength, and the 
angles of observation. If there is no dust, the diffuse component is due 
to a layer of large particles and a comparison of observations with the 
model will yield relatively small values of absorptivity (because the 
particles are large). If dust exists, the comparison will yield large 
values of absorptivity (because the particles are small). The product of 
absorptivity and particle size can be closely estimated from measuremeits 
of the amount of diffuse light. Accordingly, computations of absorptivity 
and size using measured data is one means of identifying dust presence. 
Another means may be in matching the shapes of photometric and polarimetric 
curves. 

It is interesting that the Moon possesses a two component (in size) 
optical surface. The reasons why this has occurred are not at. all clear. 
It would probably be a mistake to assume that two-size optical layers have 
occurred elsewhere, since all size combinations appear possible. In the 
case of small asteroids, there is reason to believe that dust created by' 
meteoritic impact might not be retained by the low gravity. 

Calculating Fine Dust Effects 

The p*esence of fine dust particles adhering to large particles not 
only has interesting consequences concerned with the capture of light rays. 
It also brings to mind a concept put forth in a poem by Johnathan Swift 


(1865), later paraphrased by Fielder (1920). 


"Big fleas have little fleas to 
Plague, perplex and bite *un. 

Little fleas have lesser fleas 
and so on.... ad infinitun." 

The application of this concept to airless planets greatly reduces the area 
available for specular reflection. For, if one takes the area of 

out jutting particles as 25 % ^ and this area in turn is covered with 
smaller particles, and so on, then the net area of specular reflection is 
.25 X .25 X .25 ... ad infinitun. In the model here, the effect is to 
increase diffuse light at the expense of singles, while double reflections 
remain about the same. 

It is probable that the lunar surface reflects relatively small 
amouits of singles. This would be' reflected in polarimetric data 
as a decrease of P (Max) and increase in P (Min.). Bowell, 

DollfUs and Geake (1972) have compared the optical properties of luiar 
samples with adherent dust against the polarization found for laboratory 
materials and noted the d^rease of amplitude of maximun polarization 
P(max) accompanied by an increase of P(min) anplitude in the lunar samples, 
as would be expected fhom this concept. 

In this model, fine dust is accomted for by a factor DUS which has a 
value 1.0 for total dust coverage, and 0 for none. Note that FORTRAN-style 
notation, such as "DUS" will be used for some mathematical variables in 
this article, especially those that actually appear as such in the FORTRAN 
program given in the appendix. 

The Single Reflections of the Model 

Assume that the incident sunlight is unpolarized. After each 


reflection from a surface, Fresnel's laws state there will be a component 
of light, Ej_ , polarized perpendicular to the plane of incidence, and a 
component, E//, polarized parallel to the plane of incidence. The 


reflection coefficients for the two components are taken as 

Xs2|Exp/I, = [a-b|^/la+b|^ and (1) 

Y = 2|E//f/I^‘ = |c-b|^/|c+b|^ (2) 

v^ere a = cos i 

b = (m? - sin^i)^^^ (3) 

2 2 

c r m COS i = m a (4) 

and i is the angle of incidence, equal to half the phase angle G- , and m is 
the complex index of refraction defined by m = n(1-i?(. ), and is the 
intensity of incident sunlight. 

The polarization fraction P is given by 

P = (X-Y)/(X+Y) (5) 


The singles intensity function is composed of two terms. An 
opposition term depending upon reflections from pits and affected by 
shadows, and a term due to reflections, external to pits, unaffected by 
shadows. Both terms can be diminished by the effects of fine adhering 
dust. 

For the opposition term, it will be assumed there exists an averagje 
uidth/ depth ratio, WID, which characterizes the pits of the surface, as in 
Figure 2. At phase angles G<WID, most of the surface is in shadow and the 
value of the term falls to near zero. Using an exponential to represent 
the function, the term is. 


- 10 - 


opposition intensity^ V^D(1 - DUS) exp( - G/WID)^ 


( 6 ) 


The opposition term becomes zero for singles if the insides of the pits are 
covered with dust that absorbs or converts all the incident light into 
diffuse light. 

The externally reflected term requires choosing that fraction of the 
surface which is exposed and available for specular reflection. This 
fraction, EXT, is expected, as discussed, to be in the range 0.1<EXT<0.3 
for meteorite- Impacted surfaces, but conceivably could be different for 
exotic surface materials. 

Summarizing, the singly-reflected light intensity is given by 


= IjIeXT (l-DUS)SHAD + WID 0 - DUS)(exp(G/WID)^Jj£_+JfJ (7) 

A Coronent on Reflection of Singles 

The Essence of intense reflection of light fhom the Moon and Mercury, 
at near grazing angles of observation is Prima Facia evidence that single 
reflection is .small . because Fresnel's equations predict reflectivity 
(for singles) of 40% to 100% at phase angles of 160° to 180°. This strongly 
suggests that the factor EXT (external specular area) must be^ small 
in the computer model, or that DUS (fraction of fine dust coating) is 
large or a combination of both. 

On the other hand. Mars, the Asteroids, and other outer planets and 
moons have never been observed at large i^ase angles, so the possibility of 
significant specular reflection of singles exists. The asteroids, 
especially, may be strong candidates to show some specular reflection, 
since some of them may be still yomg since formation by molecular 
condensation, and their small gravity may not retain fine dust created by 
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nieteroritic impact. The variation of the factor EXT and DUS in the 
computer model can show predictions of photopolarimetric curves for 
comparison with observations, and so it may be possible to detect such 
bodies even though observation at large phase angles is prohibited. 

A second conclusion is that models and theories of light reflection 
fhom rough surfaces which are based upon an assuned domination by single 
reflection, are fundamentally unsouid. The results of such models should 
be used with great caution. 

The Doubly-reflected Light Function 

The computer model of Wblff (1975) which made a raany-rayed simulation 
of doiisle reflections at all possible angles was convincing and yielded 
much information on the processes of light reflection, but the extreme 
complication of the computer routine made it cunbersome and not amenable to 
analysis. Using the experience of that model, a simpler analytic 
expression has been devised in v*iich four rays, one in each of four 
quadrants, are used to calculate the properties of doubly-reflected light. 

The new expression must still preserve the important properties of 
double reflection which are: : 

1) the net reflectivity depends on the product of two reflection 
coeffients; 

2) intermediate light paths into left-right quadrants are calculated 
separately from paths going into forward or backward quadrants; 

3) shadow attenuation can be applied to the forward quadrant paths, 
in accordance with the shadow geometry of a rough surface; 

4) the net polarization calculates to zero at opposition; and, 

5) the dependence on phase angle is closely preserved. 

The four representative intermediate- path rays can be chosen in many ways, 
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but one apparent choice is that ray which points directly through the 
center of each quadrant; i.e. exactly forward, exactlv backward, left and 
right. These have the required properties, but the net reflectivity is 
about 20 % less than the whole quadrant integration, so an adjustment factor 
of about 1.20 is in order. 

In the left and right quandrants, the average intermediate path is 
always at right angles to the incoming and outgoing rays , so there is very 
little change of the average angles of incidence which are about 45°. It 
is therefore sufficient to use a fixed representative ray path to calculate 
reflectivity. 

Although the angles of incid«ice can be fixed at 45°, the change of 
phase angle causes rotation and requires the use of the matrix calculation 
of Wolff (1975): 

• • ( ^ J 

These are greatly simplified because of the constant angle 45°, so 




M/a 


Cos^ (G/2) 

Sinf(G/2) 


/ 


Sin*^ (G/2) 

Cos"^(G/2) 


vdiere X(45°) and Y(45°) are reflection coefficients at 45°, I^ and are 
the initial light components (we have chosen unity), and ^ 2 ^ and are 
the final doiisle-reflected light components in the lateral (left and right) 
directions. 

Only one such calculation is needed for each phase angle instead of 
hundreds as before. 

The matrices containing sines and cosines of G/2 accouit for the 


(s)- 


Hzo\ 


:(45°) 
|Y(45°) 


M /a 
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interchange of Fhesnel components as the phase angle changes from 0° to 

1800. 

In the forward and backward quadrants, the angles of incidence of the 
average rays at the two surfaces change in direct proportion to the phase 
angle, so the representative reflectivities become. 


Forward quadrant: ^ "^2F = 1 .2 ^X(45® 

/2Y = 1.2(^Y(45°*^G/4, MR, MI)]2 

Backward quadrant: ^ X23 - 1.2 ^X(45° “ G/4, MR, MI)^2 

^28 = 1.2 (xr45®- G/4, MR, MI))2 


(9) 


(10) 


Probability Functions for Doubles 

In addition to computing the reflectivity of the two surfaces, it is 
necessary to compute the geometric probability that two reflections will 
occvr in the surface. An assumption is used here that the probability not 
used by refraction or by singles reflection, will be evenly divided between 
doiiJles and higher order reflections - as was fouid for the case of 
hemispherical pits. Maxiraun doubles probability occurs at opposition. 

There are two terras of the maximun doubles probability; one term is 
due to reflections fhora inside of the pits, and a second term is due .to 
dust~ covered area outside of the pits. If there were no dust, the outside 
area results totally in singles. These two terms are: 

GP2 = maximun doubles = 1/2(1 - EXT) + (EXT) (DUS) (11) 
geometric probability 


Equal nunber of incident light rays must go into each quadrant at 
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opposition because of symmetry. It will be assumed that this is 

approximately true at all phase angles. To support this, one observes 

that polarization of the Moon and Mars are almost independent of the angle 
of incidence. Note that phase angle independence is not assumed - quite 
the contrary. 

The probability of doubles escape from a pit depends on the phase 
angle in different ways for the different quadrants. This is because 
of the different angular paths taken by the rays as illustrated in 
Figure 2A, B, C. The forward rays are cramped for space to make two 
reflections and therefore a narrow exponential function with half the 
angular pit width is used as in Figure 2B. The backward, left and right 
quadrant rays are not constricted as the singles, so the exponential 

expression to express thePhase angle dependence has four times the width. 

Summing up, we can now write the intensity of the doubly-reflected 
light as the product of the reflectivity and the geometric probability, 
as follows: 

Left/right: Is(GP2) exp(-G/4WID) (X2 l+Y2l)/2 (12) 

Foward: l 2 p= Is(GP2) exp(-2G/WID) (X2p+Y2p)/2 (13) 

Backward: 1^^= I^(GP2) exp(-G/4WID) (X2g+Y2g)/2 (14) 

where the left and right quadrant contributions are identical. 

The total contribution of doubles is 

^2"" ^2l‘^ ^2r'*' ^2B ^ ^2F 
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II. A NEW FUNCTION TO COMPUTE DIFFUSE LIGHT 


We seek a particular mathematical description of a well-known 
property of light diffusely reflected £Y*om particles. It is commonly 
observed that colored crystals appear white if they are finely ground. 
This is because light paths through the crystals are shorter when finely 
groixid, and the total surface area is larger. Accordingly, the reflection 
opportunities, which appear white, are increased and the absorption 
opportunities, which jroduce color, are decreased. This i^enomenon is an 
example of the calculation sought. Specifically, we want to find an 
expression for the emitted diffuse light as a function of particle size and 
the optical indices of the particles. 

The advantage to be gained fhcm this approach is that by using the 
measured albedo and an estimation of the real index, a calculation can be 
made of the p-oduct; 

(particle size) x tmaginery index) 

Then, two measurenents at different wavelengths, displaying different 
absorption, will yield an unambiguous particle size. Further measurements 
confirm the result. 

This approach consists of four steps. First, the amomt of light 
which enters the particulate surface is calculated. Second, the distances 
travelled inside the particles are calculated. Third, the amount of light 
emitted into the source semisphere is obtained, and finally, a normalized 
aigular distribution function must be added to provide isotropic (assixned 
random) angular emission. 

Step 1: Light Entering the Surface. 
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The anount of light refracted into the surface particles can be 
computed by assuning that all angles of incidence are equally likely, which 
is reasonable for a random jagged surface. 

The weighted average coefficient of refraction, designated REF, is 
the integral overall angles of incidence, a: 



(16) 


which can be nunercially integrated easily in a computer frogJ"®* For 
purposes of analysis, a simple approximation can be used. 


REF 1 


- 1/2 



MI, MR)+(Y(45°, MR, MI) 


(17) 


where the integral has been approximated by the va].ue of the integral at 
the center (45®) of the integration range. Ihis is perhaps no more than 
5-1 5X in error. 

Step 2'.and 3: Light Paths Inside and Escape from the Particles. 

After a photon enters a particle, it will boince around inside by 
internal reflection, until it either escapes by striking a wall at an 
incident angle which is less than the critical angle, or, it will be 
absorbed. It follows that the length of the path traveled inside depends 
on the real index of refraction through the critical angle, on the size of 
the particle, and on the imaginary index vhich governs absorption. 

The probability of escape depends on the size of the "Cone of escape” 
whose apex angle is the critical angle and which is formed around a 
perpendicular to the wall. If all angels of striking the wall are equally 
probable, then the escape probability is the ratio of the escape-cone 
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solid-aigle to the solid-angle of a hemisphere (2tt radians) . • However, 
since an escaping ray can be emitted into either the outside hemisphere 
(really escapes), or into the planetside hemisphere (gets trapped again), 
the true probability of escape is 1/2 of that ratio. The escape cone solid 
angle will be dnoted CONE(n). For an index n=1.6, C0NE/4 tt — 125. 

The emitted ray intensity at the Nth reflection is 

% s - In cone/4tt^:i n (i8) 

By regarding N as a functional variable, with little error, the above 
expression integrates immediately to 

^N = lN-1 exp(-N C0NE/4TT) (19) 

If the subjjtance of the particle is absorbing, then the ray will al,’'- 
be diminished between reflections by absorption while travelling the paths 
of length N x D, where D is the average distance between reflections - 
about the same size as the particle diameter. This can be written 

= lo exp (-KDN) (20) 

where 1^ is the intensity of a ray vhen it first enters the particle and K 
is the absorption coefficient; K = 4 tt(mI)/X 

Combining Eqns (18, 19, 20) yields the light emitted at each 
reflection 

^^N = lo exp(-N(KD + CONE/Jfn-)) x (C0NE/4TT )L$i 
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The total «nitted light is obtained by sunming this expression over all 
reflections from 0 to oO . Again, to avoid a cunbersorae procedure and 
result, an integral can replace the sum with little error. 



o 


exp(-N(KD + CONE/4 tr)) x ( CONE/4 tr)dN 


Integrating and replacing with IgREF, yields 
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I 3 = Diffuse emitted light = IgREF/d + 4 tiKD/C0NE) 

= l5REF/(l+ (4n')^D^^MI/C0NEj (22) 

where I 5 is the illuminance of the sun at the surface of the planet . D^is 
particle diameter in wavelengths, and REF is Eqn. (16 or 17)^ already obtained. 

Equation (22) is the expression sought. It provides the intensity of 
diffuse emitted light as a function of the particle size D, and the indices 
of refraction, MR and MI which are contained in CONE, REF, and K. This 
result is somewhat surprising. One might expect an exponential attenuation 
of light to yield a sharp darkening as absorption is increased rather than 
the simple inverse relation obtained. But the sunming of many paths 
results in a less severe attenuation. 

Step 4; Angular Distribution of Diffuse Light from a Planet. 

It is evident that diffuse light from a planet travels the least 
distance through the planet's surface in reaching an observer, when it 
comes fhom the sub-earth point and is at zero piiase angle. At larger phase 
angles, the path through the surface assunes greater and greater values. 
The derivation of a function which describes the attenuation due absorption 
along the path is giv«i by Vfolff (1975, p. 1402). The result is classical 
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and termed the Lommel=Seeliger function (L-S). In terms of the phase 
angle, G, and the longitude, L, it is; 

ATTENUATION = 2 cos(G-L)/(cos(G-L) + cost) (23) 

In order to approximate the opposition effect, one can arbitrarily 
divide the light of Eqn. (24) into halves, one half of which is subject to 
opposition shadowing, thus Eqn. (24) is multiplied by 

Opposition factor = 15 + ^exp (G/WID) 

We can now reach our final conclusion by combining the above with 
the new result shown in Eqn. (22), and thus obtain a complete expression 
for diffuse light from a rough, particulate planetary surface: 

PLANET"S DIFFUSE LIGHT = Ig 2 REF cos (G-L) {h + %exp(G/WID)) 

(1 +41TKD/C0NE) (cos(G-L) + cosL) 

(24) 

The physical interpretation of Eqn. (24) is that light enters a 
particulate surface and along its path successively enters various particles, 
bouncing around inside of them, exits, enters, etc. and eventually a part 
of the light reaches the surface and goes to the observer. Two photon 
capture probabilities are involved, one is capture by a particle and this 
probability increases with phase angle, G, and depends on L; the second is 
capture by the crystalline lattice inside the crystal and this depends on 
particle size, its length, and the indices of refraction. 

Eqn. (24) is a new result. Its ultimate usefulness and accuracy in 
describing the emission of diffuse light remains to be tested by 
observation. 


Case of An Entire Planet 


If the entire planet is in the field of view of the measuring 
instrument, as must be the case when observing small objects, such as 
the unresolvable asteroids, the total signal is the integral over the 
visible, illuminated portion of the planet. A spherical body will be 
assumed. 

The differential element of area in the integral must be a projection 
of the planet's surface onto a plane perpendicular to the line to the 
observer. Reference to Figure 3 provides an element dA = (R cosL cosl dL) 
(R dl(l - cosL(l-cosl))). R is the planet radius, L is the longitude and 
1 is the latitude. The illuminated part of the planet lies between 
longitude - Tr'/Z and longitude yl'/Z-G. 


The integrand must sum over the three components of light I^, 


and I : 
3 


I(planet) = E(sun 




+ I + L dA 
2 3 J 


-jr ' % 

where the three components are given by Eqns, (7), (15) and (24). 

The integration over latitude 1 is straight-forward and results 
n " Cr 

I(planet) = R^i di ^(cosL)+cos^L(77?2 - 2)J + l^J (25) 


There appears to be no ready integral for Eqn. (25). If I = f(L), 
so in the computer model, this integration is carried out numerically. 
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Terminator of the Planet 


At the location of the terminator, a function of 6 and L, the 
intensity of 1^, and Ig must go to zero. Eqn. (23) accomplishes this 

for I^. For 1^^ and Ij, an approximation for the visible, illuminated area 
of a particle is useful : 

SHAD = 1 + 

This is used as a multiplier cif and Ig. 


adilO-cosLfi-cos^)) 



FIGURE 3* Integration over th« visible 
illuminated surface of a planet uses 
an area element projected onto a surface 
perpendicular toiho observer. 

Diffusely-emitted light has a different 
angular dependence than singles or doubles 
and has a greater dependence upon longitude. 
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F 




III. COMPUTER PLOTS OF THE NEW MODEL 

The following sets of figures use the new surface model to plot 
quantities which are observable, as well as the separate components of 
polarizatif^n and light intensity which are not observable but of interest. 

At first, the main goal of the plotting program was to produce a 
model which matched the Moon, since we know most about that surface. This 
was done by choosing suitable values of the parameters DUS, EXT, real 
index (MR), imaginary index (MI), pit width/depth (WID) and longitude. 
After those values were once determined, they were not altered during the 
production of this entire group of figure sets, except for the variation, 
one-by-one, to illustrate the effects of each parameter, shown in the 
figure sets. 

The unusual discovery was made that by changing only MI, it was 
possible to reproduce the characteristics of much of the CSM taxonomy of 
asteroids. This is discussed at length in what follows. The quantitative 
matching of computer results with observations over a broad range of 
parainoters and measurements is quite remarkable in view of the diverse 
properties concerned and the complex nature of the functions involved. 

It is concluded that the model can be regarded as quantitative for many 
purposes. 
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L REFRACTIVE INDEX (MR) VARIED 
gure lA to Figure 'IK 
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Figure 5A» Polarization at a point on a planet, 




of a point on the planet, : 
length shift of 33 percent 
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The values of P-min and the in 
match those of the C asteroids 
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The inversion at the CM Boundary Index (.10) is 
clearly seen in this figure. 
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Polarization is very large because the diffuse 
light is almost zero, 

A LARGE POIARIZATION SCALE IS USED HERE. 
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Figur e 9 A . Polarisation at a point on a planet. 
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values of the imaginary index result in almost total 
absorbtion of all ight entering particles. 
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Figure fOH » The poiariTiatxon due to the singiy-reiiectea Figure 10 I« The polarisation due to the doubly-reflected 
component of light, alone. component of li^t, alone. 
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This figure c^grees with the observation of Lyot 
that the polarization increases when near the 
terminator. The effect is small as generally 
reported. 
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Figure ;4H, The polarization due to the singly-reflected Figur e HI . The polarization due to the doubly-reflected 
component of light, alone. component of light, alone. 
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IV. THEORY OF THE POLARAMETRIC NEGATIVE BRANCH 


Besides indicating the presence of a pitted regolith, the negative branch has 
the potential of providing useful information, especially for the outer planets 

I 

which are only visible at small phase angles* It follows that there is a need 
to understand .the formation of the negative branch in Quantitative detail so 
as to make use of the numerical information available* This section describes 
in detail the light components which contribute to the negative branch,and, 
formulae are derived for the approximate position of the inversion angle and 
the polarization minimum* 

Figure ' 'L is a plot of the four coi^onents of light which contribute to the 



FIGURE 1 5 . The light components which make up the negative branch. 

The polarization .minimum depends on functions of the pitwidths only. But the 
polarization inversion angle depends on the pit functions in competition 
with the specular polarization from external surfaces. 
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There is a common misconception that the angle of minimum is half of the 
inversion angle and therefore measure of either provides the seme information. 

It is evident from the above figure that neither of these is true. The two 
angles are related to separate factors. 

The components in the figure are ‘those which contribute to polarization, and are 
defined with respect to the plane of vision (Sun-object-Earth) which is the plane 
of symmetry. Double-reflecting rays which go either side of the plane and then 
forward are termed sideways doubles . They are toe cource of negative polarization. 
Rays which reflect twice and stay in the forv^rd quadrant are termed forward- 
doubles . They contribute little to the positive polarisation because of rapid 
attenuation in getting out of the pits. Backward doubles are big contributors to 
positive polarization because their path makes it easy to get in and out of pits. 
The Location of P-min . 

In the figure near the location of P-min, the two most rapidly changing 
components are the sideways doubles (-) and the forward doubles (+) . When the 
rate of change of their sum is zero, a minimum occurs. The other two smaller 
conroonents can b e ignored, because they are both + , but since one is decreasing and 
the other increasing, even their small contribution is nullified, 
accordingly, the minimum occurs when 

where G is the phase angle, S and F are constants derived from pit width/depth. 
After carrying out the differentiation, taking logs of both sides, and solving 
for G, the result is f n P/c 


Thus the phase angle of P-min depends mainly on the geometry of the pits 
through the angular attenuation constants F for forward doubles, and S for 
sideways doubles. For example, if we use the values used in the computer program, 
S * h, F“|-, and VIED = ,20, then G(min) ^ x VJID * 9° , 
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This approxina tion isgord if G(min) is small, say less than about lU® • 
At larger angles, the approximation breaks down because of the presence of 
single (+) polarization which will make the observed angle small than the 
calculated approximation. This can be deduced from the figure, 

V.lth this understanding of P-min, one is immediately led to ask, ’Why do 
the asteroids have smaller pit width/depth ratios thm the Moon, Mars or 
Mercury ? and VJhat is the character of asteroid Ceres 1 which has an 
unusually low G(min) * 7° ? 

The Inversion Point 


The point of inversion occurs when the sum of all components is zero* 

From the figure this condition is 

/ \ “G/B _ “G/S 

singles polarization (EXT,DUST. .)+ e ■ 2 e 

It is not convienient to go further with this relation because of the 

complex dependence of the singles polarization on dust, external surface, 

M, I'd, as well as the pit geometry. However it is possible to establish 

a useful qualitative principle. First, establish that backwards doubles and 

sideways doubles have a fairly fixed ratio regardless of the type of surface. 

This is true, exactly, at opposition because both kinds of rays follow the 

same kinds of paths with average reflection angles at li5° for both. As phase 

angles become larger, the sideways-doubles average angles still remain near 

U5°, while the average angles for backv’ards doubles is h,B° - G/2 • If G is not 

too big, both types of double rays are affected almost the same. 

It follows from the above thatthe most significant effects on polarization 

near the inversion point are due to the intensity of singles relative to 

double s. Now, combine this idea with the geometric principle (discussed 

elsewhere) that the sum of pitted area (produces doubles) plus non-pitted 

area (produces single.') is a constant, and we arrive at the general qualitative 

principle that ; The Inversion Point is An Indicator of the Relative Amounts 


of Fitted and Non-pitted Area • If pits are relatively larger, the inversion 
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point moves to the right. If pits are smaller to the left. This 
principle is good enough to merit quantitative calibration using laboratory 
objects. 

Proof from Observation and Analysis 

Ample evidence is available to substantiate these conclusions. Moon 
measurements always show an unchanging G(min) if wavelength is varied, but 
G(C) changes. This is because G(min) depends on pit ratios which don't change. 

In the figures of this report where EXT or DUST are varied, G(0) always moves 
right if these two factors are increased, This is because they increase the 
pitted area. Laboratory samples show the same effect. When looking at wavelength 
dependent data, recall that the measure of size should be the wavelength itself, 
so that increasing wavelength tends to decrease non-pitted area. 
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V. THE COMPUTER MODEL COMPARED TO CSM TAXONOMY 

This sui’face nodel has provided results vhich may be closely compared to 
the CSM Taxonon^ of Eov’ell et al (15?6), Fxandnation of the mathematics involved 
show that this has occured because the new model provides quantitative 
expressions for both d iffuse light and specular light. The older model had 
no computation for diffuse light; it merely lumped in the required amount to 
obtain the desired reflectivity. 

Especially, the diffuse light equation, which is a new result, is dependent 
on particle size, wavelength and both the real and imaginery indices' 
of refraction. This res ’Its in a proper relation betvjeen those factors when 
computing diffuse light and v.’hen computing specular light* This important 
addition to the mooel and the resulting quantitative agreements w ith observation 
provide confidence that the moael describes a pitted, particulate surface 
characteristic of airless planets, in an accurate fashion, 

Tv;o kinds of Reflected Light 

The behavior of the new model becomes clear if one keeps in mind how the 
■".v.’o kinds of light, specular and diffuse , change as 'die imaginary index changes, 

Liffuse light is mainly due to light which enters particles, travels 
through them while undergoing several internal reflections, and eventually 
emerges after some absorbtion, /s the absorbtion constant (imaginary index) 
is increased, this diffuse light steadily decreases in amplitude event'^^ally 
reaching almost cero, hben the imaginary index is small, diffuse light 
dominates over the specular (singles and doubles) light. However, there 
arrives "i transition value when diffuse and specular light are about equal. 

Then, f’urther • „rease of the inaex allows specular light to dominate. 

The specular light basically depends on computations of Fresnel's equations 
the res’JuLts of which are modulated to account for the effects of the pitted 
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Figure 16* This figure was prepared for coir.parision with figures 1 and 2 of 
Powell et al (197S) which show color indices against albedo in 
a scattor-diagrao! of a few hundred asteroids. Those figures forn the basis 
for tho COM -^axcnor- of oster-cids. 

The agree'ent of this conputer nodel with the G,C, and M groups 
of asteroids is excellent. In the figure, as vou move along the heavy line 
beginning atths S-range, the imaginary index of refraction (absorbtion 
constant) steadily increases. 

The o_uantitative agreement with asteroid measurement is all the 
more remarkable in view of the fact that no adjustable parameters have been 
used for fitting. The entire computer model was created, and then matched at 
only one point: the Moon, at the point F-min - 1 % and reflectivity - 9 % • 

From that time onward, all calculations have been made with the 
same model. In view of the non-linear, complex functions for the two variables 
plotted above, agreement by coincidence seems impossible. 

b'hen making numerical comparisons, keep in mini that the computer 
model contains no variation of color due to molecular phenomena. If visible, the 
model material would appear as neutral shades of gray. 
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Figure \7, This figure vas prepared for comparison with figure 3 of 

Bowell et al (1978), which is a scatter-diagram of asteroids 
using similar coordinates. Their diagram displays the grouping of the CSM 
taxonomy. 

The agreement on the location of the S and C asteroids when 
comparing the computer model with measurements is very good, however the 
model's values of F-md.n are far too large in the M-range, or what ought 
to be the M-range, They run off scale above. 

The puzzle is complicated by the lack of knowledge of the 
real observed ir asteroids. It is not known whether, in fact, they possess 
surfaces, materials, or any characteristics ccmpareable wlththe assumed 
computer model. 

Hov:ever, the computer model itself maybe unrealistic for values 
of the imagiriary index in the metallic spec’ular range. The large computed 
values of F-min, 3 to 6 percent, are a result of the assumed presence of 
pits and particles made of glittering metallic material. Is this a correct 
description ? It might also be that the metallic material of the K asteroids 
is a malleable sort and doer not form the pits and particles assumed, 
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surface. V.Taen the inaginary index is small, Fresnel's equations are practically 

unaffected by changes in the index. As it is increased, reflectivity slowly 

increases, and there arrives a value, about MI = 0.1, when the specular light 

is about equal to the decreasing diffuse light. At larger values of the index, 

the specular light increases rapidly* However, it should be noted that both 

the parallel and perpendicular components of the Fresnel calculation increase, 

but unequally, so that the weaker asymptotically approaches the other’s value. 

close 

Thus, polarization gra::ally comes/to zero for very large values of the imaginary 
index. 

These calculations are strongly based upon the physical laws of optics. 

Even though one right quarrel with the details of diffuse reflection or the 
formation of the negative branch, the general trend of reflectivity and polarization 
are inevitably in accordance v.ith the model. Accordingly, the evidence for the 
model's interpretation of the CSH taxonomy is equally strong. Figures \ 6 and 
' 7 have been prepared using the same kind of coordinates and shoving the same 
ranges of C, £ and M asteroids as displayed by the authors of the CSM taxonomy 
(Eowell, etal, 1976) . The coincidence of these diagrams, and other numerical 
agreements, can hardly be accidental especially in view of the complex variation 
of the diffuse and specular light, and the fact that tnere are no adjustable 
parameters to create "fits" . 

The Albedo Rules 

Rules have been variously proposed and used to obtain the albedo of 
unresolveable planets using the slope of polarization, maximum polarization, 
and minimum polarization as a f’inction v’f the phase angle. Figures I g, and 
have been prepared to shov: how these rules appear from the calc’ulations of the 
computer model. In all three cases, the computer model produces a quantitative 
agreement with the empirical rules in a linear range which corresponds to the 
S asteroids and part of the C asteroids. 

These diagrams can be used to deduce reflectivity in the usual way, or 
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POLARIZATION SLOPE, % per Degree, at Inversion Angle 

Figure | . The Con^^uted Polarametric Slope versus the Computed Albedo, and some 

Comparisons ;jith Asteroid measurements. 

The computer model agrees well with the measurements of asteroids, and 
supports the contention of Zellner that the Slope-albedo rule is linear ; in the 
S-range and part of the C-range, There is a saturation of the polarimetric slope 
in the M-range, due to the transition of reflectivity dominance from the diffuse 
mechanism to specular reflection. 

The small diagrams in the lower left and the upper right show the magnitude 
and direction of corrections for conditions of diffuse light and real index 
other than '.hat assumed for this figure. 

The fact that the correction vector for point-particle diffuse light is 
parallel to the linear portion of the Slope-albedo Rule means that calibration 
error due to this cause will not occur, but subsequent use-error id.ll be concealed 
i.e. C-objects with extremely fine (d<;X) dust vd-ll be misinterpreted too far 
towards the S-range, 

There is a theoretical minimum albedo of 2.5$ for pitted, particulate bodies. 
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Figur e 'i 9 . Maxinrura and I-Iinimura Polarization From the Computer Model. 

These two curves show a theoretical version of' the corresponding rules to 
obtsin albedo from F-inax and P-min, These rules are not as reliable, it is thought, 
as the Slope-albedo rule. The reason, discussed elsewhere, is that errors in the 
P-min and F-max rules are opposite and tend to cancel in the Slope Rule. 

The computed P-max Rule computed here agrees well with observations, however 
the F-min Rule here is drastically violated in the metallic M- range by asteroids 
presumed to be of this type. 

Reasons for the P>range failure could be : the pitted, particulate model 
does not describe objects in the M-range; or, there do not exist any M-type 
objects, i.e. metallics with large absorbtion and pitted surfaces. 
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conversely, if the reflectivity canbe measured by an independent means such as 
radiometry, then the variation from the computed curve can b e used to deduce 
variations from the assumptions of the computer model. Potentially computable 
variations include the real index of refraction and the existence of very fine 
dust vhich causes diffuse light by point-particle diffraction. 

The Xaximum and Minimum Polarization Hules 

Figure ^ shov;s the computed rules. These rules have been found empiracally 
to be less accurate than a Slope-albedo Rule and the model offers a logical 
reason for this. 1-Iaod.mum polarization of the the positive branch depends on 
the amount of singly-reflected light compered to the diffuse light. The singly- 
reflected light, in turn, depends on the amount of specular area visible on the 
planet. Any variations of that area fraction must produce errors in the rule. 
Similarly, the minimum polarisation rule depends on the amount of doubly- 
reflected light ccmparec to the diffuse light, The doubl 3 ''-reflected light, in turn, 
depends on the amount of pitted, particulate area visible on the planet. Any 
variations of that area fraction proauces errors in the minimum rule. 

Thy, ■J’olarizatior. Complementarity Principle 

A very useful principle for interpreting polarimstry can be derived from 
a simple geometric concept concerning airless planets. It car. be stated, 

"A gecmeury-depenaert Increase of Positive Polarization is Always Accompanied 
hy a Decrease of Negative Polarization” . It follows from these three 
assumptions : 

1) On a planet, the area of pits plus the area of non-pits is a constant 
(equal to the total area). 

2) Pits produce double reflections and negatively-polarized light. 

3) Non-pits produce single reflections and positively-polarized light. 

Thus, if the surface structure of a planet should vary from some norm which, 
say, characterizes an albedo rule, then the singly end doubly reflected light 
vrill vary inverselj% An albedo rule dependent on the positive branch, will 
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have an error, but the albedo rule which depends on the negative branch vrill have 
the opposite error. Another example : If a geometry-dependent variation should 
occur in the study of polarization inversion angles, then a change which increases 
the positive branch will decrease the negative branch and the inversion point will 
move to the left . Oppositely, a variation decreasing positive polarization mil 
cause the inversion a igle to move to the right . 

The Slope-albedo Rule 

In a review of albedo rules, Bowell and Zellner (197ii) clearly show that 
the empirical scatter of data is least for the Polarization-slope Albedo Rule. 

The con^uter model provides a logical reason why this should be true. It results 
from the Polarization Complimentarity Principle and the fortuitous choice of the 
inversion point to measure the slope. This point lies midway between the domains 
of the positive and negative polarization branches. If geometry-dependent variations 
from the Rule norm should occur, then bj- the Principle, opposite chaz-iges occur in 
the negative aixl positive polarization '''ranches i.’hich tend to cancel each other 
in the slops. If the reader is not convinced, a little pencil and paper work 
c’rav.dng the slopes is suickly convincing. You will also notice, in pencil 
diagrams, that although the slope remains constant, the inversion angle moves 
to t he left or right , depending on the type of variation chosen, as required 
by the Principle* 

Lespite these advantages of the Slope-albedo Rule, examination of figure IQ 
shows that the . linear rule may fail completelj’’ for some dark G objects and 
metallic objects. The double branch of the computed rule suggests that some 
K or C objects could be confused with each other, if independent means is not 
used to determine on which portion of the curve an object, lies. 

A hew Reddening-Albedo Rule 

Figure 20 shows a plot of reddening vs reflectivity as computed using the 
surface model. Like other rules, it has a linear portion and a second branch. 
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REDDENING, magnitude per degree 

Figure 2.0 > Reddening is defined as the increase of color ratio with 
increasing phase angle. The data for this figure was 
obtained by measuring the slope of the computed color 
ratio curves. Reflectivity is takan at zero phase an^e. 

All data refers to curves for the integrated whole planet . 

It is interesting that in both the diffuse reflection region 
and the specular region, the rule of thumb, "The Darkest Planets ^how 
The Greatest Reddening" , is observed. However, there has never been any 
recognition of two branches to the rule, nor a failure of the rule in 
the M-range of asteroids. Reddening measurements have very large errors 
which is why it has been only a qualitative rule. Recent improvments in 
the sensitivity of photometry might change this situation# 
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It also probably is relatively free from geometry- caused error, like the 

Slope Rule, because reddening depends on the relative slopes of the diffuse 

light compared to the slope of the sum of singly and doubly-reflected light • 

If the reddening measurements are made at angles where the anplitudes of 

singly and doubly reflected light are about equal, then geometric variations 

in the surface will not produce errors. Such equality only occurs for the 

very-large values of the imaginary index and at phase angles near zero, thus 

it might, be useful for outer solar-system objects which are metallic. However, 

this is some error cancellation for all values. 

Unfortunately, the reddening rules requires great sensitivity and 

accuracy im the photometry, so that it would be useful only for bright objects. 

An Improved Polarization-minimum Rule 

The P-min ^Hile for finding albedo is very useful for objects in the outer 

solar-system where phase angles are always small, but it lacks the error-cancelling 

feature of the Slope-albedo Rule, However, an improved P-min Rule can have the 
c 

error canelling feature by using the inversion angle as a negating parameter, 

A 

For example, an adjusted P-min is calculated using 


AE JUSTED p-:aK 


fi + 

: ) P-raiiri 

k kGg J 



where is an average (calibrated) inversion angle and k is a calibration 
constant. This makes a correction for the geometry-induced errors in the 
same way as occurs naturally in the slope-albedo rule. The shift of the inversion 
angle G(C) detects the change of geometry and is used for correction. 
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COLOR RATIO (difference in magnitudes) , coir?3uter model. 

For a wavelength shift of 33 percent. 

Figur e 2 ! » This figure shows the color ratios taken from the computer plots 
displayed in figures elsewhere. This color ratio is similar to’ 
color indices such as the U-V or the r/B index of Bowell et al (1978), 
The color ratios and reflectivity numbers are taken from the 
cong:uter-generated figures at zero phase angle. 

The small diagrams at the top of the page provide corrections 
to the curves for values of parameters different than those 
assumed in the computer model. Note that the corrections are vectors. 
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Obtaining Added Information From Albedo ■t^ule Variations 


If independent means is available for measuring albedo, then the variations 
from the the values expected from polarization or reddening measurements, can be 
used to calculate other parameters of the surface. Two important such parameters 
are the real index of refraction and the amount of very- fine dust which causes 
Jiffase light by diffraction. 

Each of these two parameters causes variation of the computed reflectivity 
with a unique vector on the appropriate diagram. All portions of the diagrams 
can be calibrated vath these vectors. Tvro examples are shown on the figures. 

The ?■ and E T^'pe Asteroids 

*he R and E asteroids do not appear to fit immediately into t he series of 
objects corresponding to an ordered variation of the imaginary index. However, 
this conclusion is ver^a tentative because ihere are only 2or 3 of these objects 
identified, nevertheless, such objects probably fit into the computer surface 
model because the only deviations of the objects are with respect to color, and 
color due to molecular/atomic absorbtions is not modelled. Furthermore, the 
position of these objects in the high albedo range means that there is considerable 
diffuse light and therefore color variations are strongly present. 

Another Model Correspondence 

■‘■he computer model agrees ;vith the GSK taxoncm;' in yet another vjay. In the 
model there are various logical means to add color for S,X,E and ?. objects, but 
no Tv-ay for the dark C t-’pes. And, indeed, there is only one observed C object, 
among -the hundreds T-rhich shoTcs color. Asteroid 65 lo has a small pyroxene band 
cbsorbtion. This observation seems m;;7sterious. 
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Introduction 


VI. THE WAVELENGTH DEPENDENCE OF LIGHT 
FROM DARK, ROUGH PLANETARY SURFACE 


An obvious means of deducing some properties of the surface of airless 
planets is measurement of reflectance variation with varying wavelength. 

Such a measurement at one viewing angle will provide C'^lor information but 
cannot distinguish a flat sheet from a rocky landscape. To obtain information 
about the surface geometric structure, it is intuitively clear that the surface 
must be viewed at different angles seeking some variation as a consequence o’f 
the surface geometry. In practice, this means obtaining color data at different 
phase angles. 

For this discussion, tvo sources of color varistion are considered : Cne, 
where the absorbtion of light inside the surface material is due to v;avelength- 
dependent interactions of lifht with the atomic and molecular energy levels. Two, 
where the observed color variations are to size and g eome trie factors of a pitted, 
particulate surface such as that of the airless solar-system bodies. 

The former variations are often seen as absorbtion bands in the reflectance 
spectra. Such measurements and their analysis have been extensively carried out 
by Mike Gaffy and Tom McCord (1^76), and Clark Chapman, to name a few investigators 
in an active field. The results can lead to identification of some of the mineral 
tj-pes making up the surface. This computer model does not model this effect; instead 
it is ''neutral gray". 

The second source is not so obvious because the color variations are mainly 
due to the relative sizes of particles compared to wavelengths, Observeable color 
changes are only significant when particle sizes are of the same order as the 
vravelength, say, Several phenomena can occur. In particles which 'trap' 

light teirorarily by internal reflection, short wavelengths are absorbed more 
thffl long wavelengths, because the short wavelength path is a greater number of 
wavelengths and that is the factor upon which absorbtion depends. The result &s 
increasing diffuse light reflectance mth increasing wavelength. 
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Another phenomenon of the second type is the change of the ratio, of 
diffusely-reflected light to specular light as the viewing angle changes. 

This happens beicause diffuse light tends to emeege evenly at al.most all 
angles, i.’hereas single and double reflected light lends to emerge at an^es 
close to the incident rays, if the surface is a pitted one. Since the diffuse 
reflectivity is wavelength dependent and the speculftt* light is much less so, 
the result is a wavelength dependence on the phase angle. 

ibese latter two geometric effects are implicitly included in this 
computer noael cf a pitted particulate surface. 

Explanation of The Reddening Effect 

An unexplained phenomenon of the airless planets is the reddening of 
their reflected light with increasing phase angle, ■‘•his has been observed 
by Gehrels, Coffeer. and Owlrgs (I96I4) and Gehrels and Owings (1962) on the 
Moonj on Mars by Irvine et al (1966) and on the asteroids by V/ocley et al 
(1955). ■‘■hese observers measured reflected light using photometers and band-pass 
filters in the IR, visible, and UV wavelengths. Reddening was r eported as an increase 
of the ratio of long wavelength intensity compared to the shorter wavelength 
intensity, as a difference of magnitudes. 

Various explanations of reddening have been suggested including a) Mie 
scattering from transparent particles of a dominating critical size, b) surface 
roughness concentrated in wavelength-sized dimensions that produce suitable 
diffraction, and c) electrostatically suspended d'ust. No quantitative analysis 
has resulted. 

This computer surface model has been used to calculate reddening as a 
result of geometric effects (but not due to molecular colors). Each of the 
figure sets elsevjhere in this report has a computer plot of color ratio as 
a function of phase angle. These plots agree perfectly with the observations 
to the extend that observational errors permit (errors are large). It can be 
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. concluded tnat the geo”;etric redcening computed is sufficient to explain 
the observations of the Moon and asteroids. There is a possibility that 
molecular reddening may also increasethe effect on Mars. 

Reddening to Investigate Surface Structure 

As mentioned above, reddening measui'ements vhich take place at sucessive 
view angles contain information about the three-dimensional surface struct..H*e. 
In principle, the same is true for polarization, and polarization aifference 
measurements, for that reason, all of these plots have been provided in tnis 
report. Admittedly, I don't have any good scheme yet to use the polarization 
differences. 

There is a cor.verse iaea. If the surface structure can be ass'umed, and of 
course this ass'umpticn is implicit in all the albedo rules which rely on tne 
existence of a pitted particulate surface, then it is possible to deduce the 
reflectivity of the surface. The polarization-albedc rules are all well known, 
tut reddening has never been used in a quantitative way. However, there is a 
qualitative rule of thumb, "The Darkest Objects Always show The Greatest 
Reddening" . I believe that the improvements in sensitivity of photometry 
now make it possible to obtain sufficiently accurate data, 

A Reddenir.g-aluedo i^ule has been plotted in Figure 2.0 using the slope of 
the color ratio, versus reflectivity. Slope is taken at 30 degrees phase 
angle and reflecuivity is taken at zero phase. There are not enough accurate 
measurements existing to make a comparison with this theory meaningful, but 
it may be hoped that certain geometric sources of error will be balanced out, 
just as occurs in the polarization slope albedo rule. 

Computer Calculations 

do obtain figures (B),(C),(C) of each figure set, the following formulas 
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vero used to calculate the ordinates : 

RSDDS*^ PCBIT - 2.5 x ln(B 3 _/B 2 ) 

REDDENED PLANET - 2.5 x ln( Integra te^E^ / integrated Bg) 

POLfRIZAlTON DIFFERENCE « P^ - P, 

g 1 

vhere the subscript 1 refers to a calculation with particle size(XTL) ■ 15 
wavelengths, and subscript 2 for particle size(XTL) ■ 2C wavelengths, Ej and 
&2 are the total reflectivities. 

The corresponding wavelength shift is (20 - 15)/15 ■ 33^ or 
(20 - 15)/2C » 25^, depending on your point of view. 
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VII. THE RINGS OF SATURN 



This computer model combined with yolsrimetry by Lyot (1929), UoUfus(1979),snd 
P.Oohnson et el (1979) provides confirmation of some of the models of the 
rings which are reviewed by Morrison (chapter 12, Planetary Satellites, J. Burns 
Ed, 1977) and by Cook and Franklin (ibid, 'Chapter 19) and, yields a calculation 
of the diameter of the bodies composing the rings. 

Optical , infra-red, and microwave measurements of radiation plus radar 
returns have led to the identification of water-ice as a constituent of the 
ring bodies together xdth methane or ammonia compounds in the form of dirty 
rough snow balls. Size estimates are between 3 and 3C centimeters (Morrison) 
or 6 to ICO centimeters (Cook and Franklin), Alternatively, the balls could 
be much larger with smaller, imbedded, dense nodules of those dimensions, in 
order to account for the radar results. 

The polarimetry , measured at phase angles of 6 degrees or less, shows 
a negative branch amplitude of 0.5 to 0.6 % . The angle of minimum is ly 
defined at about 1 ± | degrees (Lyot) and 1.2 i I degrees(Johnson) . The polarization 
v;ere caused by a pitted surface, then the pits appear impossibly deep. Nevertheless 
the opposition effect has been observed with the same narrow width. The 
conclusion seems inescapable that the shadowing of opposition is respon?- 
sible for the negative branch, (The suppression of forv;ard-going double 
reflections by the shadovang is the cause of the negative branch.) 

However, it is not necessary to invoke shadowing by surface pits only. 

Inter-ball ;,hadov:ing is possible and likely, and creates no paradoxical 
geometries, so that negative polarization from it follows logically. The 
equivalent of pit/depth ratio, but for inter-ball shadowing, is the 
dimaeter/separation ratio’. Let us estimate, using the polarization data, 
that the computer variable EID Of 1 ± t degrees, or that separation/diameter 
= IL + 3 . 

It now becomes poosible to calculate the size of the dirty snowballs. The 
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reported optical depth of the rings is slightly less than 1 . The thickness 
of the rings has been measured at about 1 kiloneter. Thus a little geometry 
yields the average diameter of the balls as a little less than 


average ball diam « Tt* x RING THICK'ISSS 

I 4 . (separation/diam)^ 


28 * 10 centimeters 


This value agrees with the limitation imposed by the reported radiation 
and radar measurements mentioned above. It does exclude the very large 
snovjballs vath imbedded dense nodules. 

It is possible that there is some negative polarization due to the pitted 
surface of the balls themselves, but it is not likely. The dimensions of the 
particles and pit walls are tiny - obviously much smaller thsn the balls 
themselves. In view of the reported high albedo, it is unlikely that these 
pits and particles produce much shadowingjbecause they are translucent. Shadows 
are the sin qua non of the negative branch. If there is no shadovjing, there is 
no negative polarization. On the other hand, large balls are quite opaque, even 
though composed of particles which are translucent on a micro-scale, and 
produce first-class shadows. 

The following figure shows an estimated inter-ball polarization together 
with a small possible surface polarization which sum tothe observed data. 



Figure . Estimates of the polarization components of Saturn's 
rings which match the observed data. 
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X. THE FORTRAN PROGRAM 
(Used for Figs. 8A to 8K) 


•t • 
C _ 


■3 • 


tLt 

3. 

5. 
6 • 


-SYF PUTsICHS, ♦ISDPCH ' ■ - - 

^^OPtlS POLAPl 

REAL L* MR, PI _ 

COMPLEX CA,CE,CC,CM " " * ' 

CIMEf'jSION AB11(92,6)»AE12(92,6),AB21<92»6>,AB22(92,6>,AB1<52»6>, 

* AE2(?2,fc) , AE3<92,6),AB(92,6),AP0L(92,6),APB1(92,6),APB2(92,6}, 

* AFLND(92,BT, APL&F(92,B) ,ARED(92,6),APEEP{92,6r,AP0LyE(92i'Er7 

3 F0PMAT(1X,F6.3,15G7.2) 

A FOPMAK IX, 13, 1X,16G7.2) 

C IMTIAL CCNSTAKTE + START CO LOOPS ’ 


1 1 . 


U I C “ . 2 C 

12. 


DUSr.*- 

13. 


L=-1.T 

lA. 


f«' P = 1 . G 

15 . 


EXT=.3 

It* 


MI=.0C7 

X i • 


DC ICl hiDX = l,5 

1 ? . 


MI = .2+(f''DX-D* 

i®. 


M=NPX+1 

fi. j • 


C-P2= ( ( l.-EXT ) * 

21. 


C'-' = CMPLX(MP,-M 

2 2* 


SL'M = C .C 

2 3. 


DC 5 LR=1,9G 

2A . 


CG=(MP-D.B)/57 



CA=CMFLX(COS(D 

2 6# 


CE=CSGRT(Cw**2 

27. 


CC=CA* (Cr**2) 

2?. 


FRESY=CABS( CC- 

2 ^ • 


FF ESX = CAPS ( CA- 

3 C • 

r 

SUM=SUM+ (FPESX 

31. ■ 


REF=1 .0-SUM/9C 


iO.C) 


CB)**2/CABS(CC + CB)**'2 


C0KE=3.1A/CABS(CM)**2 
CO ICO LG=1,90 
GA=2.*(\G-1.0) 

G=GA/57.3 

PUT G VALUES INTO THE ARRAYS 


37. 

A P E C ( L' G , = 0 . t 

3 fc • 

ARECP(IvG,M)=0. 

39. 

i511(NG,l)=GA 

i~ G . 

AB12(f'.G,l) = GA 

^1. 

A521(\G,1)=GA 

^2 . 

A 5 2 2 ( N G , 1 ) = C- A 

3 . 

AB1(A‘G,1)=6A 

in. 

A E 2 < N G , 1 ) = G A 

^5 . 

A 5 3 ( \' G , 1) = G A 

Ab. 

Ac (^'G,i)=GA 

^7. 

A P 5 1 ( K G , 1) = C- A 

^ E • 

A P B 2 1 1'.' G , 1 ) = G A 


■ AF0L(NG,1)=GA 

50. 

APLMU\G,1>=GA 

51. 

AFEDPf^G,l)=GA 

52 • 

AREC(HG,1)=GA 

53 • 

APCLlvD(NG,l) = G 

5A . 

IF(G+L-1.5A) 9 


(s> E V 


HKPRODUCmiLITY OF THF 
ORIGINAL PAGI'J IS POOR 
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1 100 ISO PROPRIETARY 


=5. TO IOC 

f6. = .C^JTI^'UE 

t7. ShAD=l.+G/ (L-1.57) 

5c. ■* “ ^ C0LGr^2/Cl + C0S(Lr/C0SC0+L)) 

f=. SIHG=SIN( 0. 5*G) 

6C. CCHG=COS(0.5*G) 

61 . C FPESNEL FUhCTIOr.' RCUTINE 

62. GO TO 2 

63. 1 CA=CWFLX(COS(CG),0.0) 

eA'.'" CF = CSrvPT(C’-l**2-{STN(CG) )**2) 

65. CC = CA*(C(^**2) 

66. FPESX=CAES(CA-CB)**2/CA5S(CA+C5)**2 

67. ^ FFESY=CABECCC-C3)**2/CArS(CC+C“)**: 

63. C-C TC(10.2D»30»AC)»’JF 

69. 2 CCKTINUE 

TC. XG = .Te5A^" . 

71. ^F=l 

72. GO TO 1 

•73. lu ‘ XA5 = PFESX"‘ • 

7A. Y-i5=FcrsY 

75. C CALCULATE SINGLE REFLECTION 

■ 7&V G ®1 -TEXT* f r; -UUS) +-(1 . -EXTT* r 1 . -DUS ) *ErP f- CO/WTOTT/eTl *SHATT 

77. DG=0.5*G 

7B. NF=2 

79. - GO TO" 1 • 

30. 20 CONTINUE 

3 1 . 5X 1 = G- 1*FPESX 

Z2~, BTT=‘G^l*FPr sr~ 

33. El=0.5* (BX1+5Y1) 

?.A. C SIDEWAYS DOUBLES AND MATRIX 

S5.- — Ciy = CCOHG**21*XA5+OCIHG**21*Y45 


36. 

3'^. 




Q 


5. 
6 . 


-9 . 


' 1 


: ^ . 


6 . 

7 . 


D1Y=(SIWG**2)*XA5+(C0HG**2)*YA5 

C2X=D1X*XA5 

D2Y = DlY*rA5" ' 

X2L=(SIHG**2)*D2X+(C0HG**2)*D2Y 
Y2L=(00HC**2)*D2X+(SIHG**2)*D2Y 
C CALCULATE FCRWAPD DOUBLES 
lG=.7£5A+G/A. 

N’ F = 3 

" GC TO 1 " " ■ 

3C X2F=F=ESX**2 
Y2F = FP ESY**2 
r BACKWARD DOUBLES 
DG=.7P5A-C-/A. 
r; F = 4 

GO TO 1 • 

40 X2B=FPESX**2 
Y25=FFESY**2 
SOB = EXF(-. 13*C-/WIC) 

S:S = E>'F (-.25*G/WID) 

3X2F = GP2*EXP{-(2.0*G/WlD):i*X2F*ShAD 

EY2F=GP2*EXP(-(2.C*G/WID))*Y2F*SHAD 

B>;2L = 2.*GP2*S0S*X2L*SHAD 

BY2L=2.*GP2*S0S*Y2L*SHAD 

By 2B = C-P2*SC5*>'2B*SHAD 

5Y2B=GP2*SCB*Y2B*ShAD 

5X2=(rV2L+BX2F+BX2B)*0.3e 
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0P16AA 
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100 ISO PROPRIETARY 


D 



PY2=(=Y2L +3Y2F + BY2B)*0.38 ' - ' 

. E2=(BX2+BY2) *.5 
C FED CALCULATION 

E = .l ■ 

PLANB=. 1 
P0L=.l 

00 =6 IFEr=l«2 

XTL=n5.-^5.*aPED-l))*<l.-DUS) + l. 

E3=C0LGF*PEF*( .5+.5*EXP(-(G/WID) ) )/(l .+XTL*MI*158./C0NE) 

PEr = 2.5*AL0G(E/(Bl + B2 + B3) ) - 

B=Bl+52+B3 

C PLANETARY INTEGRATION 
SUN=C. 

RLI=-l.e7 

9 7 T" = 53M2./(1.+C0S(RLI)/C0S(RLI-^6)>)/C0LGF + P2 + B1*(1.+G/(RLI-1.57)> 

♦/EHAD 

CRC = T5*{(2./3.14)-^COS(RLI) + <0.5-2./3.14)*C0S(RLI)**2) 

SUP-SUM +DRD 

RLI=FLI-^. 04 ' " ■ ' 

IF (RLI+G-1. 57)97, 97»98 
9f. SUf = SU'-'*. G4 

REDP = 2.5*ALCG(FLANB/SUM) ‘ . - 

PLANE^SUM 

RE1=(EX1-BY3 )/(2.*E) 

RE2=(EX2-3Y2)/(2.*r) 

?0LWD=RB1+PB2-P0L 

PCL=PL1+PE2 

99 CONTINUE - 

C CO'-'PUTE FINAL RESULTS AND PUT INTO ARRAYS 
A51KNG,M) = EX1 

AB12 (NG,M ) = Evi ^ ^ - 

A 5 2 1 r : G ♦ P ) = E X 2 
A 5 2 2 ( ^' G , M ) = E Y 2 

ABl (NG,M) =E1 - ■ - - ■ ■ 

A B 2 ( N G , P ) = B 2 
AB3(NG,M) = £3 

A B ( N G , M ) = F BEPRODUeF^'i’' • OF tlffi 

A?B1(NG,P)=FB1*100.C ORIGINA' ’ ' ' 'it 

A - 52 (NG,P)=FB 2*100.0 

AF0LCNG,P)=F0L *100.0 

APLN5(NG,P)=PLANB 
APE2P(NG,P)=REDP 
A F E 0 ( ^' G , P ) = F E C 
A P G L W Cl ( ! i' G , P ) = P C L U D * 1 0 0 . 

W='ITE(6,H)riG,Bl,RED,REDP,53,B,FEl,PE2,P0L,P0LWC,PLANB,EX2L, 

*E?Y2L,“X2F,BY2R,BX2BfBY2B“ ' ^ ■ ' 

irC CONTINUE 
1 : 1 C 0 N T 1 N U 

G PLOT cpOGRAMS - - - 

CALL PLOTS ( C. , 0. , 7) 

CALL FACTCP ( 0 . ‘^ ) 

CALL PL 07(2.0, 2. 0,-3) ^17' 

A!=0L('~'1,1) = D. 

AP0L(92,1)=3C. 

CALL AX IS ( 0. , C. ,21FPHASE ANGLE - D E GR EES , - 21 , 6 . 0 , 0 . , 0 . 0 0 , 3 0 . ) 

CALL AXISiC.0,0.0,21HPCLARI2ATION RER C EN T , 2 1 , 7 . 0 , 9 0 . , - 1 0 . , 1 0 . ) 

CALL SY-'^CL ( 1 . ,6. , . 14, 3GhIPAGINARY INDEX VARIED-P RANGE ,0,30) 
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:^c;32 WGL'^F Cei6AA 


ISD PROPRIETARY 


1^9. 

17C. 

171. 
1 72 V 

172. 

I 174. 

175. 

I 17&. 
i 177. 

! I'^ST 
|l79. 


IS 4 . 
1?£. 
1ft. 
IfTT 
1 4 c . 
1 f 5 . 
i7r.'- 

191. 

192. 


CALL SYVOCL<2.0»5.Bt.l0»30HXTL = 15-20' VIAVELENGTHSTHS " *0»30) 

CALL SYK5CL(2.0»5.6».10»3CHREAL INDEX = 1.6 ♦0*30) 

CALL SYr50L(2.»5.4».l»2PHIMAG INDEX =0.2 TO 0.6 t0.t28) 

CAL L ■ STT-' BDL C 2 . 0 V 5“. 2 V. 1 0 i 2 2 HPTT'W rDTH/DEPTH"^~D'i'8 0TDT22T 
call SY^*EOL(2.0♦5.0^.10^24HDUST = 40 PERCENT *0*24) 

CALL SYVBCL(2.0f4.6».10»30HFLAT SURFACE = 30 PERCENT t0*30) 

CALL SYr50L(2.»4.6«.l»3'CHL0?';GITUDE = -74 DEGREES '■V0»30> 

DO 51D MLINE=2,6 
APCL{91,X.LINE)=-10. 

■APCLT52»FLI1VE>^'0V 

CALL Lir.E( APOLdfl) »AP0L(1»MLINE) t90tl»4tMLINE) 

CCN-^INUE 

CALL FLCT(12.0,0.0«-3) " 

A P C LW D ( 9 1 ♦ 1 ) = 0 . 

APOLWD< 92d ) = 3C. 

'C ALT.- AXTS'l 0 .‘t C . » DTETPHA SIT" A'NGLE"~DTGRETSV='2rrerV0T0 . « OTD'OT3‘OT) ’ 
CALL AXIS(0.0»0.0»21HPOLARIZATION D IF F-P C T » 2 1 » 7 . 0 » 90 . * -2 . * 2 . ) 

CALL $YM50L(l.t6.«.14»3CHIMAGINARY INDEX VARIED-M RANGE tOfSO) 

C ALL S YF B D L ( 2 . » 5 . 8 ♦ . !’» 3 0 H W A V ELZ^GT H“S H IF T ■'= " 33P C T » 0»3 0 ) ' 

DO 519 MLIND=2i6 
APOLWD (91 »MLINE )=-2. 

"AF0LWIH9 2'irL"ITJE7^-2-i 

CALL LINE(AFOLWD(l»l)»APOLWD(l»MLirvE)»90»l»4*MLINE) 

CONTINUE 

CALL PLOT( 12.0»0.0t-3)' ■; 

ARED(£4,1)=C. 

AP ED( 85»1 ) =30. 

'TALL AXIS rO . » 0 . , 21HPHASr'ANGLT''7=-TrEi;FETST'-21"»‘6VI!'rO-. f'0-TO'D»3ir."T 
CALL AXIS(C.tO.»2f;HREDD£NED POINT* DIFF-MAG » 28 ♦ 7. * 90 . * i 0 » . 05 

CALL SYMBDLd. *6. » .14,30HIMAGINARY INDEX VARIED-M RANGE *0*30) 
CALL SYMBOL (2. *5.8 ». 1 * 30HV;AVELENGTH SHIFT = B^PCT *'0»30T 

DO 517 M,LINE=2»6 
APED(a4.MLINE) = 0. 0 

AR ED( B5 *ML IN'El-i 05 

CALL LINE(ARED{1»1)»ARED(1»MLINE)*83»1«4*MLINE) 

CONTINUE 

CALL PLCT(12.C»0.t!*-3) ^ 

A!? EDP ( £4, 1 ) =0 . 

ApEDP(85*1)=50. 

CALL" AXTSrC.*0V»21HPHASE ANGLE' - D EDF EES i- 2 IT'6 .‘0 VC . V 0 i O'O »"3 0 .-) ' 

CALL AXIS(0.tC.»28HPEDDENED PLANET* DIFF-YAG * 28 * 7 . * 90 . * . 0 * . 05 
CALL SYMBOL (1. * 6 . *. 14*30HIMAGINARY INDEX VARIED-M RANGE *0*30) 

CALL SYMBCLf2.*5.f*.l*3CHWAVELENGTH SHIFT r-SSFCT" *0*30) 

DC 51£ KLINE=2*6 
APEDP(S4,MLINE. ) = .C 

' AF EBP rB'5 *7-’:LIT'IE ) OB 

CALL LINECARECPd*!)* AREDP (1 * M L I N E ) * S 3 * I ♦ 4 * ML I NE ) 

CONTINUE 

CALL PLOTC12.C*0.0*-3-) ■ 

AFLN5 ( 91* 1 ) =C, 

AFLNE(92*l)=3C. 

CALL AXIS C C . * 0 . * 2 IHPHASE AN GL E ‘ - "DEGREES * - 21 * 6 . 0 VO ." * 0 .‘0 0 * 30 .’) ' 
CALL AxIS(0.C*0.0*21HREFLECTIVITY*PLANET * 2 I * 7 . 0 * 9 0 . * 0 . 0 * . 0 5 ) 
CALL SYMBOLd.»6.*.14*30HIMAGINARY INDEX VARIED-M RANGE *0*30) 

DC 509 MLINE = 2*6 ' ' 

AFLN6(91*MLINE) = 0. 0 
A P L N 3 ( 9 2 * M L I E ) = . 0 5 
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lOD ISO PROPRIETARY 


26* 

27. 

2.‘-. 



hC . 

“ 1 • 

42. 

43. 

44. 

45. 

46. 
:47. ■ 

4f., 

4=. 

• 

«. u • 

31 . 
52. 
:53.- 
54, 
:55. 

’ c r, 

:57, 

:5P. 

s cr C - • • ■ ■ 

. w « • 

■60. 
:6i. 
•62. 
?63. 

'64, 

?65. 

■ ^ ^ 
.67. 

:f fl. 


:7C. 

'" 1 . 

:^2. 

2 73 . 

1 i ^ • 
27'=. 



- 1. 


5 09. 


511 


512 


514 


a5 


CALL LINE(APLNB<ltl)»AFLNB(l»MLINE)*90»l»4*MLINE) 

COMir.'UE 

CALL FLDT(12.0»0.0t-3> 

A53( = l»l) = 0. 

AB3(92»1)=30. 

CALL AXIS ( 0. » 0. »21HPHASE ANGLE - DEGREES » -2 1 » 6 . 0 * 0 . t 0 . 0 0 * 3 0 . ) 
CALL AXIS (C.0t0.0t21HREFLECTIVITY» DIFFUSE » 2 1 t 7. 0 ♦ 90 . ♦ 0 . 0 » . 05 ) 
CALL SV!/E0L(l.t6.».14»3CHIMAGINAPY INDEX VARIED-M RANGE «0»30) 
CO 511 KLINE=2»6 

AE3( = 1»KLINE) = 0.0 ' 

AB3(92»PLINE)=.05 

CALL LINE«AE3(1,1)»AB3(1»MLINE),90»1»4»MLINE) 

CONTINUE ' ■ • - 

CALL PLOT(12.C»O.C»-3> 

A P ( “ 1 ♦ 1 ) = 0 . 

AE(92»1)=3 0. - ■ 

CALL AXIG(C.»C.,21HPHASE ANGLE - DEG.s E ES » - 2 1 » 6 . 0 » 0 . * 0 . 0 0 ♦ 3 0 . ) 
CALL AyiS(C.r).0.0t21HREFLECTIVITY»TGTAL ♦ 2 1 » 7. 0 t 90 . t 0 . 0 * . 05 ) 
CALL SYKFOL(l.t6.».14»30HIMAGINARY INDEX VARIED-M 'RANGE tOiBO) 
DC 51? MLINE=2»6 
AE (91.MLINE) = 0 .0 

AE {92fMLINE) = .05 ' ' ' ' 

CALL LINE(AB(l,l),AB(l,MLINE)»90»l»4iMLINE) 

CONTINUE 

CALL FLOT(12.0»0.0»-3) 

APBl (rl.l )=C. 

A = E 1 ( 2 ♦ 1 ) = 3 0 . 

CALL AXISCO'. »0.»21HBHASE ANGLE OEGR EES » - 21 » 6 . 0 ■» 0 . ^-O'.tJff 3 0 . ) ~ ' 
CALL AXIS(0.0»0.C,21riPCLAPI2ATICN»SIN6LES » 2 1 » 7 . 0 » 9 0 . » 0 . 0 » 1 0 . ) 
CALL SYF5OL(l.«6.f.l4,30HIMAGINARY INDEX VARIED-M RANGE »0*30> 

DC 513 MLINE = 2»6 * - 

APB1(51»'>'LINE) = 0.C 
AP51(92.MLINE)=10. 

CALL LlNE(AFBl(l»l)»APF.l(l,MLINE)»90Vl,4,MLINEr " ' ' ’ 

C C N T I N U E 

CALL PLOT(12.C»O.C»-3) 

A » E 2 ( 9 1 , 1 ) = C . ■ 

A P P 2 ( 9 2 » 1) = 3 0 . 

CALL AXIS(G.,G.,21HFHASE ANGLE - DE GR E E S » - 2 1 ♦ 6 . 0 » 0 . » 0 . 0 0 » 3 0 . ) 
CALL AXISC0.D»0. 0»21HPGLAPI7ATTON,DOUBLES » 2 1 fT. 0 t 9 0 . » - 4 . , 4 . 0 ) 
call SyN5CL(l.»6.».14,3CHIMACINARY INDEX VARIED-M RANGE *0»3C) 
DC 514 r>^LINE = 2»6 
AF52 ( 91 .ML INE) =-4. 0 
A»E2l‘-'2,MLINE)=4. 

CALL LINE( AP£2( 1. 1 ) » AF32 ( l.MLINE) »90» 1»4»MLINE) 

CONTINUE 

call plot ( 12.0.0. Ct-3) 

A&1(91»1)=C. 

A51{=2.1)=3C. 

CALL AXIS( C. , C. .21KPHASE ANGLE - D E GP E E S . - 2 1 ♦ 6 . 0 . 0 . » 0 . 0 0 » 3 0 . ) 
call AXIS(0.0.G.C.21HREFLECTIVITY.SINGLES .21.7. 0.90. .0.0.. 004) 
CALL SYMEOLC 1 . .6. . . 14. 30HIKAGINARY INDEX VARIEO-M RANGE .0.30) 
DC 515 MLINE=2.6 
A 2 1 ( 9 1 . F L I N E ) = : . G 

An 1 (=’2.MLIi\E) =. 004 

CALL LI.‘:E(A51(1,1).AB1(1.MLINE).90.1.a.mlINE) 

CONTINUE 
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254. 
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32 ;-'OLFF 0816 A A 


C 1 


16 


C816 


1 


100 


ISD PROPRIETARY 


^fc2(91,rLlKE)ro.o 
^?2(_52 ?MLINE) = .004 

COMI WUe" 

Sfi;*’ 12. 0» 0. 0*999) 

EVD" ' • 

STOP ■ ■ ■ 
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